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Abstract. Building on recent work in SU(7V) Yang-Mills theory, we construct a 
manifestly gauge invariant exact renormalization group for QCD. A gauge invariant 
cutoff is constructed by embedding the physical gauge theory in a spontaneously broken 
SU(A'^liV) gauge theory, regularized by covariant higher derivatives. Intriguingly, the 
construction is most efficient if the number of flavours is a multiple of the number of 
colours. The formalism is illustrated with a very compact calculation of the one-loop 
(3 function, achieving a manifestly universal result and without fixing the gauge. 
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1. Introduction 

One of the hallmarks of QCD is the qualitatively very different behaviours observed in 
the high and low energy domains. In the former case, the theory exhibits asymptotic 
freedom, allowing phenomena to be described in terms of weakly interacting quarks 
and gluons. Since the coupling is small, calculations can be performed in perturbation 
theory. However, as the energy scale is lowered, so the coupling strength increases, 
ultimately causing the quarks and gluons to be bound together into hadrons. The 
failure of perturbative techniques to capture this behaviour presents a stern challenge. 

A promising approach to extracting information from the strongly coupled (or 
non-perturbative) domain of quantum field theories is the exact renormalization group 
(ERG) [1-3], the continuous version of Wilson's RG. The central feature of the ERG 
is the implementation of a momentum cutoff. A, in the theory in such a way that 
the physics at this scale — which is encoded in the Wilsonian effective action, S\ — is 
described in terms of parameters relevant to this scale. The ERG equation determines 
how evolves with A, thereby linking physics at different energy scales. Consequently, 
an ERG for QCD has the potential to provide access to the strongly coupled regime. 

In addition to providing a powerful framework for addressing a wealth of non- 
perturbative problems in a range of settings (see [4-12] for reviews), a particular 
advantage conferred by the ERG is the huge freedom inherent in its construction [13]. 
In the context of gauge theories, this freedom can be exploited to construct manifestly 
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gauge invariant ERGs [14-19] (for a comprehensive review of the alternatives, see [20]). 
Whilst being of obvious novelty value, manifest gauge invariance also provides both 
technical and conceptual benefits. From the technical standpoint, the gauge field is 
protected from field strength renormalization and the Ward identities take a particularly 
simple form since the Wilsonian effective action is built only from gauge invariant 
combinations of the covariant derivative, even at the quantum level [15]. In addition, 
the difficult technical issue of Gribov copies [21] is entirely avoided. Conceptually, a 
strong case can be made for manifest gauge invariance being the natural language to 
describe non-perturbative phenomena, not least because all conclusions drawn will be 
completely gauge independent. 

The majority of work into this scheme has, so far, focused on constructing [14- 
19], testing [16-18,22] and refining [23-25] the formalism. Most recently, however, real 
progress has been made in understanding how to compute objects of particular interest; 
specifically, the expectation values of gauge invariant operators [26]. Moreover, crucial 
steps in this procedure have a non-perturbative extension [27] . Given these developments 
we feel that it is timely to extend the framework to incorporate quarks, in anticipation 
of application to non-perturbative QCD. 

The strategy we adopt is to build directly on to the SU(A^) Yang-Mills construction, 
which we briefly describe. Recall that the implementation of a gauge invariant cutoff 
comprises two ingredients [28]. First, we apply covariant higher derivative regularization. 
However, as is well known [29], this is insufficient to completely regularize the theory, 
since certain one-loop divergences slip through. The solution we employ is to instead 
apply the covariant higher derivative regularization to a spontaneously broken SU(iV|iV) 
gauge theory, into which the physical SU(A^) gauge theory has been embedded. The 
heavy fields arising from the symmetry breaking act as Pauli-Villars (PV) fields, 
supplementing the covariant higher derivatives to furnish a complete regularization of 
the physical theory. 

The symmetry breaking is carried by a Higgs field, C. Upon acquiring a vacuum 
expectation value (ueu), this field breaks SU(A^|A^) down to its bosonic subgroup, 
SU(iV) X SU(A^) X U(l). One of these SU(iV) symmetries is identified with the symmetry 
of the physical gauge field, the other is identified with an unphysical field, A^, which 
we note comes with wrong sign action. (We can effectively ignore the U(l), as we will 
see later.) Besides A^, only A^ remains massless upon spontaneous symmetry breaking, 
all other fields picking up a mass of order A. As we send A to infinity, all effective 
interactions between A^ and A^ vanish and so the non-unitary A^ sector decouples and 
can be ignored [28]. 

The most obvious way to add quarks is to embed them in the fundamental 
representation of SU(A^|A^). However, as we will see, the structure of the SU(A^|A^) group 
in fact forbids this program. Instead, we embed N quarks into a field which transforms 
as a bifundamental under U(iV|iV). This essentially corresponds to a gauging of not 
only the physical colour symmetry but also, in an entirely unphysical way, a flavour 
symmetry. The unphysical fields which accompany the physical quarks are given a mass 
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of order the cutoff and so act as precisely the set of PV fields we require for the theory 
to be regularized. In this way, we are able to incorporate quarks in multiples of N , 
such that the elements of each set have the same mass. To obtain quarks with arbitrary 
masses, we now modify the Higgs sector to break the unphysical SU(A^) symmetry 
completely. 

Quite apart from the ERG, the inclusion of quarks in the SU(A^|A^) regularizing 
scheme is of interest in itself, since it allows the construction of a real, gauge invariant 
cutoff in QCD. In the case of pure Af = 4 super Yang-Mills, a dual picture of this was 
recently constructed [30], providing a concrete understanding of how the radial direction 
in the AdS/CFT correspondence plays the role of a gauge invariant measure of energy 
scale. It will be interesting to see whether the inclusion of quarks in the SU(A^|A^) 
regularizing scheme leads to a new way to introduce quarks in the dual picture. 

Having incorporated massive quarks into the regularization framework, it is now a 
straightforward matter to generalize our manifestly gauge invariant flow equation from 
SU(A^) Yang-Mills to QCD. To understand how we go about doing this, we review the 
structure of general flow equations. One of the key ingredients of any flow equation is 
that the partition function (and hence the physics derived from it) is invariant under 
the flow. As a consequence of this, the family of flow equations for some generic fields, 
(f, follows from [13-15, 19] 

where the functional parameterizes how the high energy modes are averaged over 

(and we have written 5'a as just S). The total derivative on the right-hand side ensures 

c 

that the partition function Z = J V(f e is invariant under the flow. 
Taking ip to represent a single scalar field, we can use [32] 

1/A 

where A is an ERG kernel and Si = 5 — 25", S being the seed action [17-19, 22, 31-33]: 
a functional which respects the same symmetries as the Wilsonian effective action, S, 
and has the same structure. Physically, the seed action can be thought of as (partially) 
parameterizing a general Kadanoff blocking in the continuum [13,19], and acts as an 
input to our flow equation. 

There is a deep relationship between the kernel and the classical, two-point vertices. 
Specifically, if we set the seed action classical, two-point vertex equal to its Wilsonian 
effective action counterpart, Sq'*''^, then we find that [32] 

V'^A = 1, (1.3) 

where A is the integrated kernel a.k.a. effective propagator: 

A = -A9aA. 

Equation (1.3) is the effective propagator relation [17] and is at the heart of the 
computational technique employed within our approach [17-19,22-26,32,33]. 



= o / ^-j/TtA' (1-2) 
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Now suppose that we consider a flow equation for some set of fields, rather than a 
single field. It is highly desirable to insist on an effective propagator relation for each 
individual field which means that, in general, the number of effective propagators — and 
hence the number of kernels — must equal the number fields. This observation holds the 
key to generalizing our flow equation for SU(iV) Yang-Mills to one appropriate for QCD. 

To construct an ERG for SU(A^) Yang-Mills, we use the template (1.1), covariantize 
the relationship (1.2) and incorporate the SU(A^|A^) regularizing structure [17,19]. By 
defining the covariantization appropriately (we will review this in section 4.1), we can 
ensure that an effective propagator exists for each of the broken phase fields. Note, 
however, that the form of the effective propagator relation is different to (1.3) in the 
gauge sector, as a consequence of the manifest gauge invariance: 



^o^^^{p) two-point classical vertex in the sector, carrying momentum p, 

and A^^;^^(p) is the associated effective propagator. Thus we see that the effective 
propagator is the inverse of the classical, two-point vertex only in the transverse space; 
equivalently, it is the inverse only up to a remainder which we call a 'gauge remainder'. 

To add quarks, we again use the template (1.1), covariantize the relationship (1.2) 
and incorporate the SU(A^|A^) regularizing structure, but this time the covariantization 
is designed appropriately for the quarks. To allow independent quark masses, we update 
the Higgs sector of the flow equation and modify the covariantizations in the gauge and 
quark sectors to ensure that there are enough independent kernels to cope with the 
breaking of the unphysical SU(A^) symmetry. 

Anticipating that the flow equation for QCD is most efficiently stated via 
its diagrammatic representation [17-19,22-25], we could jump straight to this 
diagrammatic form, using it to hide non-universal details such as the complicated 
form of the covariantization. However, before doing this, we will provide an explicit 
example of a valid covariantization, for completeness. Nevertheless, it should be 
understood that this is just one of an infinite number of choices which we can make 
but which, in practice, we never do: in actual calculations, we implicitly work with 
an infinite number of flow equations. The reason we can do this is that there exists 
a powerful diagrammatic calculus which enables us to perform general computations, 
almost entirely at the diagrammatic level [17,25,26], in a way such that nonuniversal 
details cancel out. Indeed, this calculus has been employed in pure SU(A^) Yang-Mills 
to give a diagrammatic expression for the /^-function, from which the universal answer 
(at least at one and two loops) can be directly extracted. This formula can be trivially 
adapted to QCD and we will use this to perform a very compact computation of the 
one-loop /5-function. 

The outline of this paper is as follows. In section 2, we review the regularization of 
SU(A^) Yang-Mills via SU(A^|A^) Yang-Mills. In section 3 we add the quarks, first seeing 
why we cannot embed them in the fundamental of SU(iV|iV) and then describing how 
we can instead embed them using a more elaborate scheme. We conclude this section 
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by showing how to give the quarks independent masses. In section 4 we review the 
construction of a manifestly gauge invariant flow equation for SU(A^) Yang-Mills and 
then adapt it for QCD. In section 5 we give a diagrammatic expression for the one-loop (5 
function which reproduces the universal result in the case that the quarks are massless. 
Finally, in section 6, we conclude. 

2. Regularizing SU(iV) Yang-Mills 
2.1. Embedding in Sl]{N\N) Yang-Mills 

Throughout this paper, we work in Euclidean dimension, D. We regularize SU(A^) Yang- 
Mills by embedding it in spontaneously broken SU(A^|A^) Yang-Mills, which is itself 
regularized by covariant higher derivatives [28]. The supergauge field, Afj_, is valued 
in the Lie superalgebra and, using the defining representation, can be written as a 
Hermitian supertraceless supermatrix (the supertrace of a supermatrix is defined as the 
trace of the top block diagonal element minus the trace of the bottom block diagonal 
element): 

Here, Aj^{x) = A^^t^ is the physical SU(A^) gauge field, being the SU(A^) generators 
orthonormalized to tr(rfrj') = 6°'''/2, while Aj^{x) = AI^T2 is a second, unphysical 
SU(A^) gauge field. The B fields are fermionic gauge fields which will gain a mass of 
order A from the spontaneous symmetry breaking; they play the role of gauge invariant 
PV fields, furnishing the necessary extra regularization to supplement the covariant 
higher derivatives. In order to unambiguously define contributions which are finite only 
by virtue of the PV regularization, a preregulator must be used in D = A [28]. We will 
use dimensional regularization, emphasising that this makes sense non-perturbatively, 
since it is not being used to renormalize the theory, but rather as a prescription for 
discarding surface terms in loop integrals [14, 28]. 

A^ is the gauge field for the centre of the SU(A^|A^) Lie superalgebra. Equivalently, 
one can write 

A, = A'11 + A^Ta, (2.2) 

where the T4 are a complete set of traceless and supertraceless generators normalised 
as in [28]. 

The theory is subject to the local invariance: 

6A^ = [V^,n{x)] + X^{x)l. (2.3) 

The first term, in which = — iA^, generates supergauge transformations. Note 
that the coupling, g, has been scaled out of this definition. It is worth doing this: 
since we do not gauge fix, the exact preservation of (2.3) means that none of the fields 
suffer wavefunction renormalization, even in the broken phase [17]. The second term 
in (2.3) divides out the centre of the algebra. The reason for doing this is as follows. 
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The superfield strength is jF^j, = i[V^, V^], out of which we construct the kinetic term 
~ str jF^j^. On account of str 11 = and the fact that 11 commutes with everything, it is 
apparent that has neither a kinetic term, nor any interactions. Consequently, if AP 
were to appear anywhere else in the action, it would act as a Lagrange multiplier and 
so we forbid its presence. The resulting 'no-^° shift symmetry' ensures that nothing 
depends on J\P and that J\P has no degrees of freedom. | This will prove important when 
we come to add quarks. 

The spontaneous breaking is carried by a superscalar field 

' D \ 



C 



D C 



2 



This field is Hermitian but, unlike A^, is not supertraceless and so it is valued in the 
U(A^|A^) Lie algebra. Nevertheless, the whole of C transforms homogeneously under 
local SU(A^|A^): 

5C = -i[C,a]. (2.4) 

It can be shown that, at the classical level, the spontaneous breaking scale 
(effectively the mass of B) tracks the covariant higher derivative effective cutoff scale, 
A, if C is made dimensionless (by using powers of A) and S has the minimum of its 
effective potential at: 

where Hat is the N x N identity matrix. 

In this case the classical action 5*0 also has a minimum at (2.5). At the quantum 
level this can be imposed as a constraint on S by taking (C) = cr as a renormalization 
condition. This ensures that the Wilsonian effective action does not possess any one- 
point vertices, which can be translated into a constraint on S [17,18]. In the broken 
phase, D is a super-Goldstone mode (eaten by B in the unitary gauge) whilst the C* 
are Higgs bosons and can be given a running mass of order A [14, 17, 28]. Working in 
our manifestly gauge invariant formalism, B and D gauge transform into each other, as 
we will see in section 2.2. 

In addition to the coupling, g, of the physical gauge field, the field Aj^ carries its 
own coupling, g2, (in the broken phase) which renormalizes separately [17-19,22]. It is 
often useful not to work with g2 directly but rather with 

a = gl/g\ (2.6) 

The couplings g and a are defined through their renormalization conditions: 

S[A = A\C = a] = ^ str Jd^'x (f^^,) ' + •••, (2.7) 

I It is tempting to try and remove A'^ from the algebra. However, we cannot do this directly 
since although SU(-/V|A^) is reducible it is indecomposable: is generated by (fcrmionic) gauge 
transformations. It is possible to instead modify the Lie bracket [34] but this appears only to complicate 
matters. 
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S[A = A\C = a] = ^^ str jd^x (Pl^f + ■ ■ ■ , (2.8) 

where the eUipses stand for higher dimension operators and the ignored vacuum energy. 
The field strength tensors in the and sectors, F^^ and F^^, should really be 
embedded in the top left / bottom right entries of a supermatrix, in order for the 
supertraces in (2.7) and (2.8) to make sense. We will frequently employ this minor 
abuse of notation, for convenience. 



2.2. Ward Identities 



The supergauge invariant Wilsonian effective action has an expansion in terms of 
supertraces and products of supertraces [17]: 

OO 1 

S=Y.- d''xvd''XnS^-:.^:{x,,---,Xn)stTX^^{Xi)---X^-{Xn) 
n=l •> 



+ 



E 



2' — s s 

^- m,n=0 ^"^'^^ 



d^xr ■ ■ d^Xn d^yr ■ ■ d^y^ 'S'^^'.-arifei'-- fo" - ■ ■ ,Xn]yi - ■ ■ ym) 



Xi ■■■Xn jVl ■■ - 



str (xi) ■ ■ ■ X^- {xn) str (y,) ■ ■ ■ {Vm) 



+ ... 



(2.9) 



where the X^' and Y^^ are any of the broken phase fields, with the and bj being 
Lorentz indices or null, as appropriate. The vacuum energy is ignored. We take only 
one cyclic ordering for the lists Xi ■ ■ ■ Xn, Yi - ■ ■ Ym in the sums over n, m. If any term 
is invariant under some nontrivial cyclic permutations of its arguments, then Sn {sm) is 
the order of the cyclic subgroup, otherwise s„ = 1 {sm. = !)• 
The momentum space vertices are written 

S^r.tiPu ■ ■ ■ ,Pn) {2nf 6 (±p}j = Jd^xr ■ ■ d^x^ e-^S'--^^^5;;.f:(a:i, 

where all momenta are taken to point into the vertex. We employ the shorthand 

Since we will ultimately be giving the fiow equation for QCD via its diagrammatic 
representation, it is useful at this stage to introduce diagrammatics for the action (2.9). 
The vertex coefficient functions belonging to the action (2.9) have a simple diagrammatic 
representation: 

(2.10) 

represents all vertex coefficient functions corresponding to all cyclically independent 
orderings of the set of broken phase fields, {/}, distributed over all possible supertrace 
structures. For example, 

(2.11) 
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represents both the coefficient functions S^^^^ and 5"*^^'*^^ which, from (2.9), are 
associated with the supertrace structures str C^C^ and strC^strC"^, respectively. (We 
have suppressed the momentum arguments.) Similarly, 

s: 

represents S-^^^^^^^ , S^^^^^^ and S-^^^^^'^^ . (There are no vertices which correspond to 
a trace of a single A^, since strA^ = 0.) 

The (un) broken gauge transformations follow from splitting Q into its block 
components 

and expanding out (2.3) and (2.4) (we are not interested in the no-A^ symmetry, here). 
For this purpose, it is useful to combine the fields and with the block diagonal 
components of Ao^- We denote the resultant fields by A^ and A'^, though note that 
sometimes A^ contributions can cancel out between terms. This gives the unbroken 
SU(A^) X SU(A^) X U{1) transformations [17] 

SAl = Dj^-u^ + d^nH^ sAl = Dl-uj^ + d,nH^ 

SB, = -i{B,uj^-u^B,) 6B, = ~i{By-uj^B,) 

SD = -t{DLU^ -LU^D) 6D = -i{Duj^ -uj'^D) 

and the broken fermionic gauge transformations 

5Al = -i[B,f-TB,) 5AI = -i{B,T-fB,) 

SB, = d,T - t{Alr - tAD SB, = d,f - t{Alf - fAl) 

SC^ = -i(Df-TD) SC^ = -i{DT-fD) ^ ■ ' 

SD = -i{CW - tC^) - 2iT SD = -i{C^f-fC^) + 2if, 

where D^^''^^ = d, — iAj^"^ are the covariant derivatives appropriate to the physical gauge 
field and the unphysical copy and the dot again means action by commutation. 

As noted in [17], the manifest preservation of the transformations for Aj^ and A'j^ 
in (2.12) protects these fields from field strength renormalization. The remaining fields 
are similarly protected, as follows from (2.13). 

The transformations (2.13) for B, and D, B, and D leads us to define [17]§ 

Fm = {B„D), (2.14a) 
={B,,-D), (2.14&) 

where M, N are five-indices [18, 19]. The summation convention for these indices is that 
we take each product of components to contribute with unit weight. 

§ Actually, these definitions differ from those of [17] by a sign in the fifth component. They are, 
however, consistent with [18, 19,22-26]. 
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Two Ward identities now follow from applying (2.12) and (2.13) to the action (2.9). 
The transformations (2.12) yield: 

quSZ^'SX-X- ..,P,q,r,...) = ^.-.'.f f^X- ■ ■ ,P, q + r, . . .) - S::^l::X. ..,p + q,r,.. .).(2.15) 

The effect of the transformations (2.13) are most efficiently written in the five component 
language of (2.14a) and (2.146). Introducing a five momentum 

gA/ = (gM,2), (2.16) 
allows us to write 

qNS:::^§l::.{. . . g, r, . . .) = s:::^l::.{. . . ,p, g + r, . . .) - s.r.^lzi. . . ,p + g, r, . . .), (2.17) 

where Y and X are the opposite statistics partners of the fields Y and X. (For explicit 
expressions see [18, 19].) An identical expression to (2.17) exists for when the field Fjv 
is replaced by F^. 

The Ward identities (2.15) and (2.17) can be beautifully combined using the 
diagrammatics: 

y \P V ^1 Y X Y X ,Y X Y 

g\^Xj/r = + - - -^-^ +■■■ (2-18) 

On the left-hand side, we contract a vertex with the momentum of the field which 
carries p. This field — which we will call the active field — can be either A^, A^, Fr or 
Fr. In the first two cases, the triangle > represents Pp whereas, in the latter two cases, 
it represents pr = {pp, 2). (Given that we often sum over all possible fields, we can take 
the Feynman rule for > in the C-sector to be null.) On the right-hand side of (2.18), 
we push the contracted momentum forward onto the field which directly follows the 
active field, in the counterclockwise sense, and pull back (with a minus sign) onto the 
field which directly precedes the active field. Since our diagrammatics is permutation 
symmetric, the struck field — which we will call the target field — can be either X, Y or 
any of the un-drawn fields, as represented by the ellipsis. 

Allowing the active field to strike another field necessarily involves a partial 
specification of the supertrace structure: it must be the case that the struck field either 
directly followed or preceded the active field. In turn, this means that the Feynman rule 
for particular choices of the active and target fields can be zero. For example, as trivially 
follows by multiplying together supermatrices, an F can follow, but never precede an 
Aj^, and so the pull back of an Aj^ onto an F should be assigned a value of zero. The 
momentum routing follows in an obvious manner: for example, in the first diagram on 
the right-hand side, momenta q + p and r now flow into the vertex. In the case that the 
active field is fermionic, the field pushed forward / pulled back onto is transformed into 
its opposite statistic partner, as above. 

The half arrow which terminates the pushed forward / pulled back active field is 
of no significance and can go on either side of the active field line. It is necessary to 
keep the active field line — even though the active field is no longer part of the vertex — 
in order that we can unambiguously deduce flavour changes and momentum routing, 
without reference to the parent diagram. 
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We illustrate (2.18) by considering contracting l> into the Wilsonian the effective 
action two-point vertex: 

^>(sy = @Jr - . (2.19) 

Given that > is null in the C* sector, the fields decorating the two-point vertex on the 
right-hand side can be either both A^s or both fermionic. In the former case, (2.19) 
reads: 

P,Sff{p) = Sf{0)-Sf{0) = 

where we note that S^' is in fact zero by itself, as follows by both Lorentz invariance 
and gauge invariance. In the latter case, (2.19) reads: 



6 



N5, 



where we have used (2.14a) and have discarded contributions which go like S^'{0). 
However, the S""'(0) must vanish. This follows from demanding that the minimum of 
the superhiggs potential is not shifted by quantum corrections [17]. Therefore, 



SJ- =0. (2.20) 

2. 3. Taylor Expansion of Vertices 

For the formalism to be properly defined, it must be the case that all vertices are Taylor 
expandable to all orders in momenta [14-16]. Consider a vertex which is part of a 
complete diagram, decorated by some set of internal fields and by a single external 
(or A"^), which we denote by a wiggly line. The diagrammatic representation for the 
zeroth order expansion in the momentum of the external field is all that is required for 
this paper [18, 19]: 



M lO 
X <)\ Y 



- "^^-^ +■■■; (2.21) 



note the similarity to (2.18). 

The interpretation of the diagrammatics is as follows. In the first diagram on the 
right-hand side, the vertex is differentiated with respect to the momentum carried by 
the field X, whilst holding the momentum of the preceding field fixed (we assume for 
the time being that both X and the preceding field carry non-zero momentum). Of 
course, using our current diagrammatic notation, this latter field can be any of those 
which decorate the vertex, and so we sum over all possibilities. Thus, each cyclically 
ordered push forward like term has a partner, cyclically ordered pull back like term, 
such that the pair can be interpreted as 



- df, 



Vertex, (2.22) 



where r and s are momenta entering the vertex. In the case that r = —s, we can and 
will drop either the push forward like term or pull back like term, since the combination 
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can be expressed as (9^; we interpret the diagrammatic notation appropriately. If any of 
the fields decorating the vertex carry zero momentum (besides the explicitly drawn A*), 
then they are transparent to this entire procedure. Thus, they are never differentiated 
and, if they precede a field which is, we must look to the first field carrying non-zero 
momentum to figure out which of the vertex's momenta is held constant. 



3. Adding Quarks 

3.1. Massless Quarks 

The simplest way to try to incorporate quarks into the setup is to embed them into the 
fundamental representation of SU(A^|A^): 

•=(:)• 

where \1/ transforms under SU(A^|A^), ■?/; is a physical quark field and ip is an unphysical, 
bosonic spinor (here and henceforth, we suppress spinor indices). || Immediately, we can 
see that this embedding is inconsistent with (2.3): the supergauge invariant quark term, 

9^ 

does not satisfy no-J^ symmetry. 

If, however, the fields \& transform as R, for some representation, R, then we 
can construct a no-A^ invariant representation simply because \l/ has zero 'charge' under 
no-A^. Thus, the strategy we pursue is to embed the quarks into fields which transform 
as a bifundamental of SU(iV|iV). To achieve this, we first embed the up-like quarks 
[up, charm, top (suitably generalized for 7^ 3)] and down-like quarks (down, strange, 
bottom) into two tensor fields {'ipu)] and {'ipd)], where the superscript indices carry an 
SU(A^) colour symmetry and the subscript indices carry an (unphysical, gauged) SU(A^) 
fiavour symmetry. In turn, ip^ and ipd are now embedded into fields and which 
are valued in complexified U(A^|A^): 

Notice that and are not Hermitian. Consequently, g (p) are not related to the 
physical fields ipu ii^d) and, since they will be seen to come with wrong sign action, 
should be interpreted as unphysical degrees of freedom. These fields, together with 0^, 
0^, (f^ and <y9^ (the components of which are bosonic spinors), will be given a mass of 
order the cutoff. Under gauge transformations, and "^d transform homogeneously: 

5^„ = -^[*„,fi], 5*, = -2[*rf,fi]. (3.2) 

II If this scheme were to work, we would also have to introduce further unphysical fields to provide 
sufficient PV regularization. 
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The SU(A^|A^) invariant quark kinetic term that we include in the Lagrangian is 

just 

-^[^ti^uf ■'^u + ^ti^df -^d), (3.3) 

where the minus sign compensates the sign buried in the supertrace, ensuring that 
the physical quark terms come with the correct overall sign. Notice that we have not 
included any covariant higher derivatives; it is straightforward to demonstrate that the 
supergroup structure, alone, is sufficient to provide the necessary regularization in the 
quark sector by repeating the analysis of [28], but this time including the fields 

To show the types of terms that we must include in the action to give the unphysical 
fields a mass of order the cutoff but leave the quarks massless (for the time being), it is 
useful to construct the following projectors: 

1/1, , / lljv \ 1 , / \ ,^ 

With a slight abuse of notation, we can write = cr+\E'ucr+, ip"^ = cr_\E'uCr_, 
g = cr_^ucr+ and tpu = ct+^uCT-. We can lift these projectors to the symmetric phase 
by defining 

,±^i(]l±C). (3.5) 

Thus, to give a mass to e.g. g in the broken phase, all we need to do is add to the 
Lagrangian the term: 

-^Astr (^„c^„^„c:+) . 

Upon spontaneous symmetry breaking, this reduces to —A/g'^tigg (plus interaction 
terms). 

Thus, in the broken phase, the only massless fields we are left with are and A^. 

In the pure gauge case, we know from [28] that decouples from A^: integrating out 
the heavy fields, the lowest dimension gauge invariant effective interaction left between 
A^ and A'^ is the square of the two field strengths (according to standard perturbative 
power counting with g in the usual place) 

A-^tr(F^)2tr(F2)2, 

which is clearly irrelevant. Adding the quarks, however, we see immediately from (3.3) 
that these fields can combine with the unphysical gauge field to form a term of mass 
dimension D, 

To remove this term, and thus ensure that the unphysical fields decouple as the 
regularization scale is sent to infinity, we must tune a (see (2.6)) to zero at the end 
of a generic calculation. In fact, this tuning is perfectly natural from a non-perturbative 
perspective: the theory carried by A"^ is, as a consequence of its wrong sign action, not 
asymptotically free but instead trivial. Moreover, we can remove the need to perform 
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any such tuning by modifying the Higgs sector to completely break the unphysical 
SU(A^); indeed, we will do precisely this when we adapt the formalism to give the 
quarks independent masses. 

We conclude this section by discussing the additional (un)broken invariances which 
arise from the inclusion of the quarks. The up-like quarks supplement (2.12) with 

5ipu = «(^V« - 'ipuUj'^) 



6g = iiu'^Q — Quj^] 
5ip^ = iuj^ ■ ip^ 



(3.6) 



(3.7) 



and (2.13) with 

6ipu = i{rLp'^ - v?V) 
5q = iijif^ — if^f) 

5<^l = i{TQ - tpuf) 

S(f^ = i{f^u - Qt) 

(similarly for the down-like quarks). Notice that the quark field ifju is not protected 
from field strength renormalization. However, the transformations (3.7) do enforce that 
all components of have the same field strength renormalization (likewise ^d)- The 
unbroken transformation for ipu given by (3.6) confirms our interpretation that the 
physical colour symmetry is carried by , whereas the unphysical flavour symmetry is 
carried by A^. 

3.2. Massive Quarks 

If we only needed to give all the up-like quarks one mass and all the down-like quarks 
one mass, then we could simply add a mass term to the Lagrangian, using the fields 
that we have already: 

— \ [m„ str + nid str ^a'^+'^d'i- 

where we have used (3.5). Of course, to give all the quarks different masses, we will 
have to break the unphysical, gauged flavour symmetry. To do this, we introduce two 
new dimensionless superscalars, Cu and Cd which, like C, lie in the adjoint of U(A^|A^) 
and transform homogeneously: 

6Cu = -i[Cu,^], 6Cd = -i[Cd,^]. 
We choose the vevs of C„ and Cd to be 

where (T„ = diag(l, mc/mu, mt/m„) and, given the unitary matrix, U, WadU = ad = 
diag(l, ms/m^, mft/mrf) (with an obvious generalization for arbitrary values of N). The 
vev of Cu breaks the unphysical SU(A^) down to U(l)^~^ We could choose the vev of 
Cd to be diagonal. However, in this case we would be forced to tune the couplings of the 
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residual U(l)s (which we note are, on account of their wrong sign action, asymptotically 
free and not trivial) to zero. Consequently, we might as well choose the vev of Cd such 
that in combination with the vev of Cu the unphysical SU(A^) is completely broken. 
Thus, we implicitly assume U to be such that amongst the generators broken by {Cd) 
are those which are not broken by {Cu)- The vevs of Cu and Cd (unlike that of C) are not 
protected from quantum corrections, which of course corresponds to renormalization of 
the quark masses. The broken Ward identities (whose modification due to the breaking 
of the unphysical SU(A^) we will discuss shortly) protect the components of Cu,d from 
field strength renormalization, in the broken phase. 

With the introduction of Cu and Cd-, there is no requirement to retain C. However, 
the following exposition is made simpler if we keep C and so we do so, noting that 
such considerations are anyway irrelevant from the point of view of the diagrammatic 
form of the flow equation. We will not give an explicit realization of the symmetry 
breaking potential V{C,Cu,Cd) which yields (2.5) and (3.8), since we are free to work 
with any potential which satisfies the following requirements. First, (in unitarity gauge) 
all Goldstone bosons are eaten by the various components of which acquire mass 
{i.e., the potential must not possess any accidental symmetries: the largest continuous 
symmetry group is just SU(A^|A^)). Secondly, the remaining (Higgs) components of C, 
Cu and Cd are given a mass of order the cutoff. 

For non-degenerate masses, it is useful (and always possible) to construct the set 
of projectors which live in the bottom right block of a supermatrix: 

Pi = diag(0^,l,0,0,...) 
P2 = diag (0^,0, 1,0,...) 

; (3.9) 
To see this, simply note that e.g. 

(1 - mt/mu){mc/Tnu - mt/niu) 
We can lift the Pi to the symmetric phase by introducing the non-degenerate (running) 
parameters Oj and defining: 

= n ^^i^^. (3.10) 

In the broken phase (which recall that, by construction, we are actually always in) we 
identify the Oj with the elements of o"u. Note that the Vj gauge transform homogeneously. 
The quarks' mass term can be taken to be 

^ [m^jStr (Cn^-^u-J+^n*?-) + m^str {Cd^,-^ d^+'^ d<i~) \ , (3.11) 

g2 L J 

where it is understood from now on that we have rotated the down-like quark fields 
to the mass basis (this is exactly analogous to the introduction of the CKM matrix 
in the standard model). The remaining components of '^u,d are given masses of order 
the cutoff. Note that we have included in (3.11) purely for convenience, to ensure 
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that the masses of the components of the fields 0^ and ip^ (see (3.1)) do not pick up 
contributions from the quarks' mass matrices. 

Neglecting the covariant higher derivative regularization, the kinetic terms for C„ 
and Cd take the form 

(V^ ■ C^f + (V^ ■ Cdf . 

Notice that this term provides differing contributions to the masses of the various 
components of B. Specifically, B decomposes into columns, with each column receiving 
a different mass. This is precisely what we would expect from the unbroken gauge 
transformations, as we now discuss. 

An immediate effect of breaking the unphysical SU(A^) is that the relation- 
ships (2.12) and (2.13), (3.6) and (3.7) (which we supplement by those appropriate 
for Cu^) decompose. The only relationships which are completely unaffected are those 
involving just oji i.e. the unbroken relationships for A^, C^, 9?^ and cj)^. 

The relationships involving just a'2 are completely broken. This means that the 
independent components of each of the bottom right block fields are no longer related 
by unbroken gauge transformations and so can be expected to propagate separately. 

In the fermionic sectors, the previously unbroken gauge transformations involve 
both uji and UJ2 e.g. 6B^ = —i{B^uj'^ — uj^B^). Now, however, the UJ2 part is completely 
broken. In matrix language, the surviving unbroken transformation involving uji mixes 
up elements of each column with elements of the same column. Consequently, upon the 
breaking of the unphysical SU(A^), B^ decomposes into 'flavours', -B^^, corresponding 
to the columns, with unbroken transformation law 

{5B^,r =i[uj'y^{B^,y . 

This is precisely what we want: in the quark sector, colour remains a good symmetry 
and the unphysical, gauged flavour symmetry is completely broken. 

With the above decomposition of many of our fields, we must adapt our expansion of 
the action in terms of fields (2.9) by appropriately expanding the set of fields represented 
by X and Y . To maintain the supermatrix structure, we should ensure that the new 
fields are still embedded in supermatrices. For example, the field B^a should be in the 
appropriate column of the top-right block of a supermatrix, with all other elements set 
to zero. Equivalently, we can project the field B^a out by using the Pi of (3.9). 

At first sight, the breaking of the unphysical SU(A^) considerably complicates 
the Ward identities. However, we can anticipate from [23-25] that we can and 
should hide these complications in the diagrammatics. This will be made particularly 
straightforward if we now sum over the flavours of the target fields in (2.18). This helps 
for the following reason. Consider (2.15) where the target fields are fermionic. We know 
that the fermionic fields decompose by column and so the right-hand side of (2.15) will 
contain a sum of terms such that, if the unphysical SU(A^) were restored, these terms 
could be combined back into block supermatrix components. 
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We conclude this section by giving the renormahzation conditions for the quarks: 



1 /" 

S=— Jd^x tr [t/j^ (i + an) + i^d (i + ^d) ^Pd 



+ ■ ■ 



9 

where the eUipsis represents all other operators contributing to the effective action. 

4. A flow Equation for QCD 
4.1. Review of S\J{N) Yang-Mills 

4.1.1. Setup We begin by describing the flow equation used for pure Yang-Mills [19], 
the basic form of which is: 

-A9A5 = ao[5,S,]-ai[S,], (4.1) 

where = g'^S—2S {S, we recall, being the seed action). "P On the right hand side of the 
flow equation is the bilinear functional, aQ[S,'Eg], which generates classical corrections 
and the functional which generates quantum corrections. These terms are given 

by 

where the A represent the ERG kernels and the notation {A} denotes their 
covariantization [14,15]. 

The natural definitions of functional derivatives of SU(A^|A^) matrices are used [16, 
17,28]: 



_S_^ f 6/6C^ -6/6D \ 
6C~\6/6D -6/6Cy' 



(4.4) 



and from (2.2) [17,28]: 



' 2r.-^ + ^-4-. (4.5) 



5A^ ^5AAf. 2iV5^o 

The wonderful simplicity of (4.2) and (4.3) arises from the realization that the fine 
detail of the flow equation (which, as we will see, does not affect universal quantities 
anyway) can be buried in the definition of the covariantization. Nonetheless, for the 
purpose of transparently generalizing to QCD, we now discuss the covariantization in 
some detail. The primary ingredient is the supercovariantization [17] of the kernel W, 
{W}_A.- This is defined according to 

u{W}^v= / Wf,,...f,„^^,...^^{xi,...,Xn]yi,...,ym]x,y) 

str [u{x)Af,,{xi) ■ ■ ■Af,„{xn)v{y)Au,{yi) ■ ■ ■ A,„(ym)] , (4.6) 

V Notice that wc use Sg instead of the Si of (1.2), on aceount of g being sealed out of the eovariant 
derivative. 
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where u and v are supermatrix representations transforming homogeneously as in (2.4) 
and where, without loss of generality, we may insist that {W\ satisfies u {W\v = 
V {W}j^u. For simplicity's sake, we have chosen (4.6) to contain only a single supertrace. 
(In the diagrammatic form of the flow equation, such details make no difference.) The 
m = 77, = term is just the original kernel, i.e. 

WX--x,y) = W,y. (4.7) 

The requirement that (4.6) is supergauge invariant enforces a set of Ward identities 
on the vertices W^^...^^^^^...^^ which we describe later. The no-A^ symmetry is obeyed 
by requiring the coincident line identities [15]. These identities are equivalent to the 
requirement that the gauge fields all act by commutation [16], ensuring that the no-^*^ 
part of (2.3) is satisfied. A consequence of the coincident line identities, which also 
trivially follows from the representation of (4.6) in terms of commutators, is that if 
v{y) = llg{y) for all y, i.e. is in the scalar representation of the gauge group, then the 
covariantization collapses to 

u{W}^v = {stTu)-W-g, (4.8) 

where we define 



f-W-g 



x,y 



f{x)W.,yg{y)= / f{x)W{-d'/A')g{x) 



which holds for any momentum space kernel iy(p^/A^) and functions of spacetime /, 
g, using 



xy 



Wi-d'/A') 6{x -y)= ly(pVA^) e 



ip-{x-y) 



At this point, it is instructive to recall the demonstration of the SU(A^|A^) 
invariance of the flow equation, assuming that the covariantizations {A} are just of 
the type (4.6) [17]. 

Under (2.4), the C functional derivative transforms homogeneously: 




5C' 



and thus by (4.6), the corresponding terms in (4.2) and (4.3) are invariant, 
functional derivative, however, transforms as [17]: 



d 



6 



5A, 



.ft 



il 

H tr 

2A^ 



6 



SA, 



(4.9) 
The A 

(4.10) 



The correction is there because (4.5) is traceless, which in turn is a consequence of 
the supertracelessness of (2.1). The fact that 6/6A does not transform homogeneously 
means that supergauge invariance is destroyed unless the correction term vanishes for 
other reasons. 

Here, no-^'' symmetry comes to the rescue. Using the invariance of (4.6) for 
homogeneously transforming u and v, and the invariance of S and S, we have by (4.10) 
and (4.8), that the A term in (4.2) transforms to 

SS 



SA 



SAu 



2N 



tr 



SAn 



A-^^-str^ 
SAa 



(^^S,),(4.11) 
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where S ^T^g stands for the same term with S and Eg interchanged. But by (4.5) and 
no-^'' symmetry, 

str T-ir- = T-TK = 

(similarly for S), and thus the tree level terms are invariant under (2.3) and (2.4). 
Likewise, the quantum terms in (4.3) are invariant and this completes the proof that, 
for covariantizations of the form (4.6), the flow equation is both supergauge and no-A'^ 
invariant. 

As it stands, the covariantization (4.6) is not general enough for our purposes: we 
require the broken phase fields to come with their own kernels. The first part of the 
solution to this [17] is to define a new covariantization 

u{W}^,v = u{W}^v - \[C,u] {WMC.v]. (4.12) 

The C commutator terms are introduced to allow a difference between A and B kernels, 
and C and D kernels, in the broken phase. They do this because at the level of two-point 
flow equations C is replaced by a in (4.12), and a (anti) commutes with the (fermionic) 
bosonic elements of the algebra. Thus, extracting the broken phase two-point, classical 
flow equations from (4.2), we find that the A^ kernels are given by A-^-^, the C* kernels 
by A^^, but the B kernel is A-^^ + A-^ and the D kernel is A^^ + A^f [17]. The B and 
D kernels can be combined [of. (2.14a) and (2.146)): 

^MN\P) - y Q -/\DD j ■ 

As must be the case, the extra term in (4.12) is consistent with both supergauge 
and Yio-A^ invariance. If u and v transform homogeneously, then so do [C, u] and [C, v\. 
Correction terms proportional to the identity (c/. (4.10)) are killed by the commutator 
structure; this structure also ensures no-^'' invariance. Notice that (4.8) holds for the 
extended covariantization. 

We are still not quite done: it is convenient to generalize the covariantization yet 
further. In particular, since the physical coupling, g, and the unphysical coupling, (^2, 
renormalize differently [17, 19], it is useful to furnish A^ and A^ with different kernels 
(there is no need to do this for and [19]). To construct a term which does this, 
we cannot use a commutator, since [a, 5/5^4*] = 0. The simplest solution is to define 
in (4.2) and (4.3) 

u{A-^]v ^ u {A-^\,v + u {A-^^],V{v) + V{u) {A^],v, (4.13) 

where 

8iVP(A:) = {C,X}strC-2CstrCX (4.14) 

V{X) has the following properties which ensure both no- .4° invariance and supergauge 
invariance: 

strP(X) = 0, (4.15a) 
P(]l) = 0. (4.156) 
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In the broken phase, we see that 
5 \ 15 



V 



5A„ 



1 



where we have noted from (4.5) that strcr5/(5^^ = and the elhpsis includes terms with 
additional fields. Thus, a minus sign is introduced in the sector, compared to the A} 
sector, and we find that A^'^' = A^^ + A^^, A^'^' = A^^ - A^^ [19]. 

We conclude this section by commenting on the Ward identities satisfied by the 
vertices of the covariantized kernels. Throughout this section, we have used the fact 
that the u and v in e.g. (4.13) are functional derivatives with respect to, say, and 
to label the kernels viz. A^^^^. The diagrammatic form of the Ward identities (2.18) 
holds for the vertices of the kernels, also, so long as two of the target fields are identified 
with the ends of the kernels i.e. with Z^ and Z^ [18, 19]. 



4.1.2. Diagrammatics As mentioned in the introduction, the most useful representa- 
tion of the fiow equation is a diagrammatic one, which is shown in figure 1 [18, 19, 23, 25]. 



{/} 



® 



{/} 



Figure 1. The diagrammatic form of the flow equation. 



The term on the left-hand side generates the flow of all cyclically independent 
Wilsonian effective action vertex coefficient functions which correspond to the set of 
broken phase fields {/}. 

The objects on the right-hand side of figure 1 have two different types of component. 
The lobes represent vertices of action functionals. The object attaching to the various 
lobes, * is the sum over vertices of the covariantized ERG kernels [15, 17] and, like 
the action vertices, can be decorated by fields belonging to {/}. The fields of the action 
vertex (vertices) to which the vertices of the kernels attach act as labels for the ERG 
kernels. We henceforth loosely refer to both individual and summed over vertices of the 
kernels simply as a kernel. The dumbbell-like term corresponds to the classical term, 
ao, whereas the padlock-like diagram corresponds to the quantum term, ai.^ The rule 
for decorating the classical and quantum terms is simple: the set of fields, {/}, are 
distributed in all independent ways between the component objects of each diagram. 

Embedded within the diagrammatic rules is a prescription for evaluating the group 
theory factors. Suppose that we wish to focus on the fiow of a particular vertex coefficient 
function which, necessarily, has a unique supertrace structure. For example, we might 

+ There is an additional, improperly regularized term generated by the flow equation which has been 
removed by a suitable constraint on the covariantization [14, 17, 19, 25]. 
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be interested in just the S^^*-^^ component of (2.11). On the right-hand side of the flow 
equation, we must focus on the components of each diagram with precisely the same 
supertrace structure as the left-hand side, noting that the kernel, like the vertices, has 
multi-supertrace contributions. In this more explicit diagrammatic picture, the kernel 
is to be considered a double sided object (for more details see [18, 19]). Thus, whilst the 
dumbbell like term of figure 1 has at least one associated supertrace, the next diagram 
has at least two, on account of the loop (this is strictly true only in the case that kernel 
attaches to fields on the same supertrace). If a closed circuit formed by a kernel is 
devoid of fields then it contributes a group theory factor, depending on the flavours of 
the fields to which the kernel forming the loop attaches. This is most easily appreciated 
by noting that stra^ = ±A^ (see (3.4)). In the counterclockwise sense, a cr+ can always 
be inserted for free after an A^, or F, whereas a cr_ can always be inserted for free 
after an A^, or F. 

The above prescription for evaluating the group theory factors receives 1/N 
corrections in the A^ and A^ sectors, as a consequence of the SU(A^) completeness 
relation [15]. If a kernel attaches to an A^ or A^, it comprises a direct attachment and 
an indirect attachment. In the former case, one supertrace associated with some vertex 
coefficient function is 'broken open' by an end of a kernel: the fields on this supertrace 
and the single supertrace component of the kernel are on the same circuit. In the 
latter case, the kernel does not break anything open and so the two sides of the kernel 
pinch together at the end associated with the indirect attachment. This is illustrated 
in figure 2; for more detail, see [18, 19,26]. 



direct 



1 



Figure 2. The corrections to the group theory factors. 



We can thus consider the diagram on the left-hand side as having been unpackaged, 
to give the terms on the right-hand side. The dotted lines in the diagrams with indirect 
attachments serve to remind us where the loose end of the kernel attaches in the parent 
diagram. 

4-2. Adding Quarks 

The game now is easy. We add to the flow equation classical and quantum terms for the 
fields and modify the superhiggs sector. For all fields, we ensure that there is 

sufficient freedom in the covariantization to allow enough different kernels for the broken 
phase fields. 

Though the details will ultimately be hidden in the diagrammatic form of the fiow 
equation, we will give examples of choices we can make for the covariantization. Just 
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as in the pure Yang-Mills case, the kernels of the propagating fields may turn out to 
be linear combinations of those which appear in the fiow equation. Though we will not 
explicitly perform this change of basis, it is instructive to see how we can, if we so desire, 
construct the covariantization so as to make this procedure as easy as possible. 

As in the pure Yang-Mills case, the starting point for constructing the 
covariantization is (4.6). To this we now add additional terms which refiect the complete 
breaking of the unphysical SU(iV) gauge symmetry. Knowing that decomposes into 
fiavours, it makes sense to add to the covariantization a term of the form 

1 ^ 

-7E([^.JC,n]]{A/^}jP„[C,t;]]^ 

where we have used (3.10). (As usual, the overall factor is merely a matter of 
convention.) The presence of the C commutators is purely for convenience. In the 
broken phase, they project on to the block off-diagonal components of u and f , which 
ensures that, at the two-point level, the above term does not interfere with the fiow of 
of the components of the field A^. This makes it easier to extract the kernels of the 
propagating fields in terms of the kernels in the fiow equation. 

The gauge field has A^^ — 1 independent components. Since we can use the same 
kernel for a field and its Hermitian conjugate, we require a total of A^(A^ -|- l)/2 — 1 
kernels. So, for the components of A^, we add to the covariantization a term of the 
form: 

Af(Af+l)/2-l 



where 



V'^{X) = {Y,XZ, + Z,Xr,)strF,Z, - r,Z,strF,XZ, - Z,y}str 

and the Yj and Zj contain linear combinations of the symmetric phase versions of the 
projectors defined by (3.10) in their bottom right block (all other elements being zero). 
There are clearly many different choices we can take for Yj and Zj. An example for 
N > 2 would be to set Yj = Zj and choose the first A^(A^ — l)/2 Yj to be such that, in 
the the broken phase, they reduce to the ^6*2 independent combinations of projectors 
of the form Yj = + Pi^k- Then we can take the remaining Yj to be any A^ — 1 of 
the A^ Pi. Notice that V satisfies the conditions (4.15a) and (4.156) and thus does not 
spoil either the supergauge or no-A'^ invariance of the fiow equation. The appearance 
of V (see (4.14)) in (4.16) is again for convenience ensuring that, at the two-point level, 
the kernels of (4.16) do not appear in the fiow of the components of Bfj_. 
Finally, then, a suitable choice for the covariantization of A-^-^ is: 

1 ^ 

uiA-'-^jv ^ u {/^^^\v - - E [C, u]] {^f^}AV,. [C, t;]]; 



N(N+l)/2~l 
3=1 
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The modifications to the covariantization in the superhiggs sector are almost 
identical; the only real difference is that, since the superscalars are not supertraceless, 
there are more propagating degrees of freedom in than in the sector, and so we must 
introduce additional kernels to take account of this. 

The inclusion of quarks follows a similar pattern. The contributions of the up- 
like quarks to the flow equation are the standard ones for spinor fields [33]*, with the 
contribution to the classical term given by 

and the contribution to quantum term given by 

where we have suppressed spinor indices and functional derivatives with respect to \& 
are defined as for any other unconstrained [i.e. not supertraceless) superfield (see (4.4)). 
The covariantization is chosen to be: 

u{k''-^-}v = «{A*"-^"}.^ + E[^.>^-M] {A;"*"}.[^.,^+(^)] + • • • (4.17) 

The first term on the right-hand side is the usual contribution involving just the 
supercovariantization. In the second term we have introduced, for convenience, the 
objects 

^±(X) = 1([C,[C,X]]±2[C,X]]) 

defined so that, in the broken phase, they reduce to a^Xa^ and a^Xa^ (plus interaction 
terms). Consequently, at the two-point level, the second term on the right-hand side 
of (4.17) contributes only to the flow of the physical quarks. The ellipsis represents 
additional terms which provide kernels for the rest of the propagating fields embedded 



4-3. Diagrammatics for QCD 

4-3.1. The Exact Flow Equation The beauty of the diagrammatic form of the flow 
equation given in figure 1 is that it can be directly generalized from SU(A^) Yang-Mills 
to QCD: all we need to do is to extend the set of broken phase fields which can decorate 
the diagrams both as internal and external fields. A consequence of this is that the 
prescription for extracting the group theory factors receives further corrections, which 
follow from inserting the appropriate projectors to go from derivatives with respect to 
supermatrix block fields to derivatives with respect to their appropriate components. 

In fact, it is useful to employ a slightly different flow equation, in which we work 
directly with renormalized fields [31-33] (recall that the components of all fields bar 

* Though, compared to QED, there is no need to expUcitly take care of the anticommuting nature of 
the quark fields, since this is automatically taken care of by the embedding into the supergroup. 
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and are protected from field strengtli renormalization) . This flow equation is 
given in figure 3. It is not the result of scaling the wavef unction renormalizations 
out of the version of (4.1) appropriate to QCD but is a perfectly valid flow equation, 
nonetheless [31,32]. This is a manifestation of the tremendous freedom we have in 
constructing flow equations, encapsulated in (1.1). 



{/} 



{/} 



Figure 3. The diagrammatic form of the flow equation. 



The term I]xe{/} 7*''^'' explicitly takes account of the anomalous dimensions of the 
fields which suffer field strength renormalization. The field x belongs to the set of fields 
{/} and the notation 7^'^-' just stands for the anomalous dimension of the field x (which 
is zero for all but the components of ^u,d)- 

4-3.2. Perturhative Diagmmmatics In the perturbative domain, we have the following 
weak coupling expansions [15, 17-19]. The Wilsonian effective action is given by 

00 1 

S = Y. a'^'^'^S, = -So + Si + ---, (4.18) 

where 5*0 is the classical effective action and the S'j>o the ith-loop corrections. The seed 
action has a similar expansion: 

00 

5 = E/'^„ (4.19) 

and the (3 function and anomalous dimensions are defined as usual: 

00 

P ^AdAg = Y.g''^'P„ (4-20) 



i=l 

00 

7 



= AdA In Z^''^ = J2 9^'li''\ (4.21) 
1=1 

where Z^^'^ is the field strength renormalization of the field of species x- 

We also introduce /3 functions for the dimensionless mass parameters, frf = m^/A: 

00 

P^ = AdAW =Y.g^'Pi. (4.22) 

1=1 

Defining Sj = Si — 2Si, the weak coupling flow equations follow from 
substituting (4.18)-(4.22) into the flow equation, as shown in figure 4. 
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. i{f} 
n 



E 

r=l 



2K-i)A, + i:3;^ 




f} 
\ 

J 



{/} 



Figure 4. The weak coupling flow equations. 



The symbol • means — A9A|m:„,mc,...- We will see shortly why the notation for the 
ERG kernels, * . includes this symbol. A vertex whose argument is an unadorned 
letter, say n, represents S'„. We define = n — r and n± = ± 1. The bar notation of 
the dumbbell term is defined as follows: 

O-^iSn-r-, Sr\ = ao[Sn-r, Sr] — ao[Sn-r, Sr] — ao[Sn-r, Sr]- 

We illustrate the use of the fiow equation by considering the fiow of all two-point, 
classical vertices. This is done by setting n = in figure 4 and specializing {/} to 
contain two fields, as shown in figure 5. We note that we can and do choose all such 
vertices to be single supertrace terms [18, 19]. 




(4.23) 



Figure 5. Flow of all possible two-point, classical vertices. 



Following [14-19,32], we use the freedom inherent in S by choosing the two-point, 
classical seed action vertices equal to the corresponding Wilsonian effective action 
vertices. Equation (4.23) now simplifies. Rearranging, integrating with respect to A 
and choosing the appropriate integration constants [18,19], we arrive at the following 
relationship between the integrated ERG kernels and the two-point, classical vertices. 

We have attached the integrated ERG kernel, denoted by a solid line, to an arbitrary 
structure since it only ever appears as an internal line. The field labelled by M can be 
any of the broken phase fields. The object ► = » is a gauge remainder (c/. (1.4)). 
The gauge remainder components are non-null only in the sectors corresponding to 
(components of) and Fr and, in these sectors, > and > are related as a consequence 
of gauge invariance, as we will see shortly. Note that, in the case that a full gauge 
remainder bites a vertex, as opposed to just a >, we can replace the half arrows on the 
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right-hand side of (2.18) (which we recall just indicate to former presence of a >) with 
a > [18,19]. 

We have been able to construct the effective propagator for each and every 
independent classical, two-point vertex because we ensured that, for each such vertex, 
there exists an independent (integrated) kernel. 

From the effective propagator relation and (2.20) follows a series of diagrammatic 
identities. In QCD, as opposed to pure Yang-Mills, the renormalization of (C)^^ — 

' 1 2 

equivalently the renormalization of the quark masses — means that one-point C„'^ 
vertices exist beyond tree level, spoiling (2.20) at the loop level. The first of the 
diagrammatic identities is, then, the classical part of (2.20): 

= 0. (4.25) 

From the effective propagator relation and (4.25), two further diagrammatic 
identities follow. First, consider attaching an effective propagator to the right-hand field 
in (4.25) and applying the effective propagator before > has acted. Diagrammatically, 
this gives 



0) = = > - »> , 

which implies the following diagrammatic identity: 

» =1. (4.26) 
The effective propagator relation, together with (4.26), implies that 

'o) > = > _ >» =0. 



In other words, the (non-zero) structure > kills a classical, two-point vertex. But, 

by (4.25), this suggests that the structure > must be equal, up to some factor, to 

<l. Hence, 

> = < , (4.27) 

where the dot-dash line represents the pseudo effective propagators of [18, 19]. 

In practice, pseudo effective propagators only ever appear in a very specific way [25], 
which we now describe. Consider a three-point vertex attached to two arbitrary 
structures, A and B, by two effective propagators. The third field on the vertex is 
taken to be an A-^ carrying momentum p and we suppose that we now Taylor expand 
the vertex to zeroth order in p. Using (2.21) we have: 
'A 

^A 

- - ^ , (4.28) 

where the arrow on the momentum derivative symbol indicates in which sense the 
momentum derivative acts. This is unnecessary in the parent diagram on the left-hand 
side of (4.28), since it is obvious that the momentum derivative has 'pushed forward' and 
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so corresponds to a derivative with respect to the momentum flowing into the vertex, 
from the structure B. However, once the vertex has been removed via the effective 
propagator relation, this information is lost, unless explicitly indicated. 

Allowing the active gauge remainders in (4.28) to act (according to (2.18)) leads us 
to define 

/: ^ A 

= - - ^ , (4.29) 

'■■ ^ i 

which we will reqire shortly. Notice that the second term of the right-hand side contains 
a pseudo effective propagator (differentiated with respect to to the momentum flowing 
into its bottom end). 



5. The One-loop f3 function 

As an illustration of the formalism, we would like to reproduce a standard result, namely 
the one-loop (3 function for massless QCD (we take the massless case since, in our mass- 
dependent scheme, massive quarks will spoil universality, rendering any comparison 
with other methods of limited use). Fortunately, due to the developments in [24,25], 
this calculation is extremely easy. In [25] a diagrammatic expression has been derived 
for the n-loop /3 function in SU(A^) Yang-Mills from which the universal answer (at least 
at one and two loops) can be directly extracted. The key to deriving this expression is 
the effective propagator relation, which we have ensured holds for QCD. Indeed, at the 
one loop level, the pure SU(A^) Yang-Mills expression is exactly the same as in QCD, 
modulo the changes to the Feynman rules, and it is given in figure 6. (Beyond one loop, 
the expression of [25] is only slightly modified.) We define □^i/(p) = p^^/^i/ — P^iPu and 
take wiggly lines to denote physical gauge fields (with Lorentz indices suppressed). Note 
that in D = 4 only the second, third and final diagrams contribute, so the expression 
for (3i is really very simple. 



4/?in^.(p) + 0(p^) = - 




Figure 6. Diagrammatic expression for f3i. 
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From comparison with the Yang-Mills (YM) expression [19,22] and the QED 
expression [33], it is immediately clear that, as in conventional approaches (though 
with the number of flavours set equal to twice the number of colours), 

^QCD ^ ^YM ^ 2N^. 

Setting the quark masses to zero yields 

.QCD ^ /II _ 4 

(47r)2 U 3 



6. Conclusion 



We have constructed a manifestly gauge invariant ERG for QCD and have used it to 
compute the one-loop (3 function for SU(A^) Yang-Mills coupled to 2N quarks. In the 
massless limit, we recovered the universal result. 

The formalism is a direct extension of the one developed for SU(A^) Yang-Mills 
in [19]. The incorporation of the quarks comprised three steps. First, the quarks 
had to be added in a way that respected the SU(A^|A^) regularization scheme. The 
symmetry associated with the centre of this algebra in fact prevented simply embedding 
the quarks into the fundamental representation of SU(A^|A^). Instead, we first embedded 
sets of quarks into N x N matrices, whose rows (columns) were labelled by colour 
(flavour). Each of these matrices was then embedded into a separate supermatrix valued 
in complexified U(A^|A^). In this way, we were able to include multiples of quarks, 
with each set of having degenerate masses. 

The second step was to give each of the quarks independent masses, and this 
required that we broke an unphysical, gauged SU(iV) flavour symmetry which is carried 
by one of the fields belonging to the SU(A^|A^) regularizing structure. To do this, we 
introduced Higgs fields for each of the sets of quarks, whose vevs are essentially 
mass matrices. The introduction of non-degenerate masses lifts the restriction that the 
number of flavours must be a multiple of the number of colours, since we are at liberty 
to remove quarks from the spectrum by tuning their masses to infinity. Clearly, though, 
the construction is most efficient when, suggestively, the number of flavours is a multiple 
of the number of colours. 

The third and final step was to adapt the ERG equation. This involved not just 
including additional terms for the additional fields but also defining the covariantization 
of the ERG kernels appropriately. The key point is that we require the number of 
independent kernels to be equal to the number of independent, propagating fields. 
Having given algebraic examples of the covariantizations needed, we wrote down the 
full QCD fiow equation in its diagrammatic form. Besides explicitly including the 
anomalous dimensions of the renormalizing fields, this expression has exactly the same 
form as that used for pure SU(A^) Yang-Mills. Indeed, this similarity allowed us to 
directly write down a diagrammatic expression for the one-loop (3 function in QCD. By 
setting the quark masses to zero we were able to recover the standard, universal result. 
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The development of a manifestly gauge invariant ERG for QCD is timely. In SU(A^) 
Yang-Mills, methods exist for computing the expectation values of gauge invariant 
operators without fixing the gauge [26] and this work can now be directly generalized 
to include quarks. We find this particularly exciting in view of the fact that important 
progress is being made in understanding the structure of non-perturbative contributions 
both to these expectation values and the [3 function [27]. 
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